NPS  ARCHIVE 
1966 
MARUCCI,  T. 


BOILING  HEAT  TRANSFER  ^aITH  LOW 
QUALITY  STEAM 
BY 
THOMAS  FRANK  MARUCCI,    LT.    USCG 

SUPERVISOR     -     ARTHUR  E.    BERGLES 

MAY  21,   1965 


'A-^S;^ 


Library 

U.  S.  Naval  Po«tHradiiate  Sctiopy 

Monterey,  California 


BOILING  HEAT  TRANSFER  VJITH  LOW  QUALITY  STEAM 


by 


THOMAS  FRANK  MARUCGI,  Lieutenant,  United  States  Coast'  Guard 
B.S/,  United  States  Coast  Guard  Academy 
^  (1959) 


suB^a:TTED  in  partial  fulfillment 

OF  THE  requirements  FOR  THE  DEGREE  OF 

master  of  science  in  mechanical  engineering 

and  the  professional  degree 

naval  engineer 

at  the 
MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGI 
May,  1965 


Signature  of  Author, 


Department  of  Naval  Architecture  and 
Marine  Engineering,   May  21,  1965  « 


Certified  by 


Thesis  Supervisor 


Accepted  by 


Chairman,  Departmental  Committee 
on  Graduate  Students 


f\r 


>c 


?^ 


(^[^■■.^'J^CC 


faibrary 

U.  S.  N*val  Poatgradeate  School 

M..,„.......     -alifornia 

ABSTRACT 


Boiling  Heat  Transfer  with  Low  Stearo  Quality 
Thomas  Frank  Marucci 
Submitted  to  the  Department  of   Naval  Architecture  and  Marine  Engineer- 
ing on  Kay  21,  196^,  in  partial  fullfillment  of  the  requirements  for 
the  Master  of  Science  Degree  in  Mechanical  Engineering  and  the  Profes- 
sional Degree,  Naval  Engineer, 

Heat  transfer  for  water  and  steam  mixtures  in  a  horizontal  stain- 
loss  steel  tube  under  non-boiling,  surface-boiling,  and  fully-developed 
boiling  conditions  was  investigated. 

Local  heat  transfer  coefficients  are  presented  for  a  limited  range 
of  flow  parameters.  Inside  diameter  of  0.2ii25  inch  and  a  range  of  lii.i; 
to  78.5  L/D  are  considered.  Exit  pressure  varied  from  60  to  about  72 
psia  for  mass  velocities  from  ^.9  x  10^,  to  2.38  x  10^  Ibs/hr  ft^.  The 
range  of  steam  quality  covered  was  0  to  UU^,  The  majority  of  this  data 
was  taken  at  inlet  temperatures  of  about  290°F.  Heat  fluxes  ranged  from 
2  X  10^  to  1.8  X  10°  Btu/hr  ft^  and  were  limited  by  burnout  conditions. 

Heat  transfer  results  for  non-boiling  and  subcooled  nucleate  boil- 
ing are  in  agreement  with  previous  investigations.   Bulk  boiling  data 
is  presented  at  the  two  lower  flox'j  rates.  This  data  shows  that  as  quali- 
ty is  increased,  the  local  heat  transfer  coefficient  also  increases.  This 
is  indicated  on  the  boiling  curve  as  a  shift  toward  lower  wall  superheat 
as  quality  increases.  Similar  results  are  indicated  to  a  lesser  degree 
with  subcooled  liquids  at  low  values  of  subcooling.  This  appears  to  be 
due  to  the  relatively  large  volume  fraction  of  vapor  created  in  the  nu- 
cleate boiling  process. 

Experimental  heat  transfer  coefficients  obtained  are  considered  to 
be  in  fair  agreement  ■with  predicted  values. 


Thesis  Supervisor;  Arthur  E.  Bergles 

Title:  Assistant  Professor  of  Mechanical  Engineering 
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NQMENCLATUHE 


A 

^ 

surface  area 

^c 

tm 

test  section  cross  section  area 

S 

- 

specific  heat 

D 

M 

diameter 

E 

- 

voltage  across  the  test  section 

h 

- 

shunt  voltage 

G 

- 

mass  velocity 

ht 

>■ 

enthalpy- 

h 

■1 

heat  transfer  coefficient 

I 

- 

current 

k 

B 

thermal  conductivity 

L 

- 

heated  length 

P 

- 

pressure 

q 

K 

heat  transfer  rate 

q/A 

- 

heat  flux 

(q/A)i 

- 

heat  flux  at  the  inception  of  boiling 

(q/A)ex. 

S 

critical  heat  flux 

(qA)Fc 

IE 

forced-convection  heat  flux 

R 

« 

electrical  resistance 

T  «  temperature 

T,^  ■  inside  wall  temperature 

T^^  "  outside  wall  temperature 

V  »  velocity 

w  «  mass  rate  of  flow 

x'  "  position  along  the  test  section  measured  from  the 
inlet  to  the  heated  section 

X  «  stesm  quality 

AT  "  temperature  drop  through  the  tube  wall 

A  Tgj^t  "  wall  -  minus  -  saturation  temperature 

A^sub  "  saturation  -  minus  -  bulk  temperature 

/\  ■  dynamic  viscosity 

(T  ■  surface  tension 

p  «  density 

J>  »  resistivity 

X  "  latent  heat  of  vaporization 

Dimensionless  Groups 

Nu  "  Nusselt  number       ■    hD/k 

Pr  »  Prandtl  number       -     C^^M/k 

Re  -  Reynolds  number      «     GD//^ 
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Subscripts 

b  ■  bulk  or  mixed-mean  conditions 

BO  ■  burnout  condition 

i  "  inlet  bulk  condition 

exit  ■  test  section  exit  bulk  condition 

L  ■  liquid  condition 

sat  "  saturation  condition 

t  ■  tube 

V  ■  vapor  condition 

w  "  condition  at  heat  transfer  surface 


INTRODUCTION 

Tne  current  development  of  high  performance  heat-transfer 
equipment  has  stimulated  a  great  amount  of  research  in  boiling 
heat  transfer.  Past  studies  of  boiling  heat  transfer  have  been 
a  factor  in  the  development  of  nuclear  reactors  and  rocket  motors. 
Vaporization  of  the  coolant  is  frequently  desired,  such  as  in  the 
boiling-water  reactor  and  bootstrap  rocket  engine „ 

There  has  been  a  vast  amount  of  work  in  the  field  on  forced 
convection-boiling  heat  transfer.  Most  of  this  has  been  concen- 
trated on  the  prediction  of  burnout  under  given  conditions  of 
fluid,  surface,  and  flow  parameters.  These  studies  have  produced 
numerous  correlations;   the  best  of  which  can  predict  critical  heat 
flux  to  about  plus  or  minus  J>Q%   of  experimental  data. 

Within  the  past  ten  jears  there  has  been  an  increase  in  the 
number  of  investigations  designed  to  study  heat  transfer  rates  en- 
countered in  boiling  heat  transfer.  The  results  of  these  studies 
have  been  presented  in  two  basic  forms,  ■'^irst,  in  the  form  of  the 
boiling  curve  v^ere  heat  flux  is  plotted  versus  wall  superheat 
(T^-Tg^^)„  This  has  usually  been  utilized  where  the  bulk  fluid  is 
subcooled,  and  indicates  that  heat  flux  can  be  represented  as: 

q/A  -  C  {T„  -  T^^^)"  (1) 

Other  studies  considered  bulk  boiling  and  present  results  that  show 
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the  local  heat  transfer  coefficient  (h)  increasing  as  steam  quality 
increases  beyond  a  certain  range.   This  range  has  been  generally 
agreed  upon  as  beginning  at  about  1'^%   quality.  As  quality  varies 
below  this  value  h  is  constant. 

Generally,  when  liquid  flows  through  a  heated  tube,  subcooled 
nucleate  boiling  can  take  place  before  the  bulk  of  the  fluid  reaches 
its  saturation  point,  if  the  impressed  heat  flux  is  high  enough.  The 
eqxiilibrium  thermodynamic  quality  would  be  zero  in  this  subcooled 
region  of  flow,  but  the  void  fraction  could  be  substantial,  depending 
of  course  on  the  fluid  properties  and  heat  transfer  rate.  As  flow 
proceeds  down  the  tube  the  bulk  of  the  fluid  becomes  saturated,  net 
yapor  generation  commences.  In  this  portion  of  the  flow,  the  vapor 
generation  can  take  place  by  additional  nucleate  boiling  and/or  by 
evaporation  at  the  liquid-vapor  interface.  If  the  tube  is  long  enough, 
it  is  possible  to  reach  the  point  of  100^  vapor. 

This  study  v;as  conducted  to  extend  the  fully-developed  boiling 
curve  i  n  the  region  of  bulk  boiling  and  to  determine  what  effect  of 
variation  in  h  with  quality  has  on  the  boiling  curve.  It  was  also 
pointed  out  by  Professor  Bergles  that  non  equilibrium  voids  created 
in  subcooled  region  should  have  an  effect  on  h. 

The  equipment  used  limits  the  study  to  low  pressure,  low  flow 
rates  and  moderate  heat  fluxes.  Fluid  conditions  included  low  sub- 
cooling  and  low  steam  quality.  This  study  is  particularly  applicable 
in  the  development  of  evaporators  for  saline  water  conversion. 
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PRQCl-DURE 
Operation  and  Data  Taking 

The  experimental  apparatus  (Figure  1)  used  in  this  study  was 
one  already  in  operation  at  the  M,I,T„  Heat  Transfer  Laboratory* 
Only  the  test  sections  had  to  be  fabricated  to  suit  this  study. 
Generally,  the  test  sections  had  to  be  instrumented  to  provide 
means  of  determining  local  pressure  and  the  temperature  at  the 
heat  transfer  surface.   To  obtain  these  measurements,  pressure 
taps  were  located  along  the  tube  and  thermocouples  were  placed 
on  the  outside  wall  of  the  test  section.  To  insure  that  the  out- 
aide  wall  was  adiabatic,  two  guard  heaters  were  used  and  main- 
tained at  a  temperature  equal  to  the  tube  wall  temperature,  A 
brief  discussion  of  the  experimental  apparatus  and  a  detailed 
discription  of  the  construction  of  the  test  section  is  presented 
in  Appendix  A,  A  diagram  of  the  heat  transfer  test  section  is 
shown  in  Figure  2.  A  preliminary  investigation  of  critical  heat 
flux  was  conducted  to  be  able  to  predict  the  critical  heat  flxix 
for  the  range  of  parameters  to  be  studied.  This  was  done  pri- 
marily to  avoid  burnout  of  the  elaborately  instrumented  heat 
transfer  test  sections.  The  test  sections  used  to  determine 
critical  heat  flux  were  not  instrumented  to  obtain  wall  temp- 
eratures and  were  constructed  with  only  one  pressure  tap  located 
at  the  test  section  exito 

Prior  to  installation,  the  test  section  was  given  a  final 
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internal  cleaning  vdth  acetone.  After  installation,  the  entire  test 
section  was  wrapped  with  insulation  to  minimize  external  heat  losses. 
The  loop  was  then  filled  with  distilled  water  and  all  lines  vented  to 
remove  air.   The  degassing  operation  was  then  initiated  by  heating  the 
water  in  the  degassing  tank.  With  the  water  boiling,  the  loop  water 
was  heated  and  elowly  sprayed  into  the  top  of  the  tank.  Results  from 
Reference  (1)  indicated  that  20  minutes  of  such  degassing  reduced  the 
dissolved-air  content  to  approximately  1  c.c. /liter.  Two  phase  flow 
can  influence  system  stability  and  recent  studies  have  indicated  that 
system  instability  can  cause  premature  burnout.  Air  trapped  in  the 
system  would  provide  a  certain  amount  of  compressibility.  Reference  (2) 
discusses  system  instabilities  in  more  detail.  After  degassing,  the 
heat  exchanger  cooling  water  was  turned  on  to  bring  the  system  to  nor- 
mal temperature.  The  loop  was  then  ready  for  operation,  circulating 
distilled,  deionized,  and  degassed  water. 

Critical  heat  flux  daoa  was  obtained  by  increasing  power  gradu- 
ally wi'iile  maintaining  constant  flow  rate,  exit  pressure,  and  inlet 
temperature.  Critical  heat  flux  was  evaluated  at  the  last  recorded 
conditions  prior  to  physical  burnout  of  the  test  section.  To  vary 
exit  quality  conditions,  inlet  temperature  was  varied  in  subsequent 
runs. 

The  general  procedure  for  taking  heat  transfer  data  was  similar 
to  that  used  for  critical  heat  flux  data.  The  major  difference  was  that 
at  each  power  level,  vjhile  other  conditions  were  maintained,  the  guard- 
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heater  power  \^s   adjusted  to  the  balance  point,  i.e.,  where  the  average 
guard-heater  temperature  (two- thermocouples)  was  the  same  as  the  aver- 
age outer  tube  wall  temperature  (three  thermocouples).  These  tempera- 
tures could  not  be  kept  exactly  equal,  but  were  maintained  to  within 
plus  or  minus  ^°F   of  each  other.  At  equilibrium,  wall  temperatures, 
fluid  temperatures,  pressures,  pressure  differences  snd  test  section 
voltage  and  current  were  recorded.   The  heat  flux  was  then  increased 
to  the  next  level  and  the  process  repeated. 

After  the  highest  heat  flux  was  completed,  the  power  was  reduced 
and  the  inlet  temperature  changed.  Readings  were  then  taken  for  a  set 
of  heat  fluxes  at  a  new  inlet  temperature. 

The  maximum  attainable  heat  flux  for  any  given  run  v;as  determined 
by  burnout  except  in  one  instance  due  to  voltage  limitations  of  the 
generators.  The  burnout  heat  flux  was  estimated  from  data  presented 
in  Reference  (3)  and  from  data  taken  in  this  study.  It  was  desired 
to  approach  but  not  reach  burnout,  -temperature  and  pressure  levels 
were  largely  determined  by  the  system.  Beginning  at  the  lowest  inlet 
temperature  attainable,  the  inlet  temperature  was  increased  to  about 
20°F  below  the  saturation  temperature  for  the  corresponding  pressure. 
The  exit  pressure  was  maintained  at  60  psia  except  at  final  runs  where 
the  minimum  exit  pressure  obtainable  for  the  desired  flow  rate  exceeded 
60  psia.  This  was  attributed  to  the  fact  that  choked  flow  conditions 
were  present  due  to  the  large  volume  of  vapor  content  in  the  test  section. 

Initial  data  taken  with  T.S.-l  and  2  was  at  relatively  low  values 
of  heat  flux  and  at  low  inlet  temperatures.   This  was  to  obtain  non- 
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boiling  heat  transfer  data  which  is  well  correlated  and  hence  serves 
as  a  check  on  instrumentation.  By   increasing  the  heat  flux  slowly  it 
was  possible  to  obtain  data  on  the  inception  of  nucleate  boiling  and  to 
be  able  to  trace  the  boiling  curve  in  the  transition  region  between  non- 
boiling  and  fully-developed  subcooled  nucleate  boiling. 

At  the  higher  heat  fluxes  with  intermediate  inlet  temperatures  tho 
fully-developed  subcooled  nucleate  boiling  curve  could  be  well  estab- 
lished. 

The  final  data  was  then  taken  with  T.S.-2  at  inlet  temperatures 
near  saturation.  The  inlet  temperature  was  maintained  below  saturation 
at  all  times  since  the  equipment  did  not  allow  determination  of  inlet 
conditions  when  net  vapor  existed  at  the  inlet.  This  procedure  pro- 
vided data  for  bulk  boiling  over  the  quality  range  from  0  to  2U^.  To 
extend  the  data  to  higher  steam  qualities  it  was  necessary  to  fabri- 
cate T,S.-3o  T.S.-3  was  constructed  similar  to  T.S.-l  and  2  except 
that  it  was  twice  as  long.  The  location  of  pressure  taps  and  wall 
thermocouples  for  each  test  section  is  given  in  Appendix  A.  Due  to 
inconsistent  results  as  discussed  later,  T,S,-3  was  modified  by  moving 
the  thermocouples  at  the  upstream  locations  to  positions  corresponding 
to  L/D's  at  the  downstream  position  of  T,S.-2.  One  series  of  runs  was 
conducted  with  this  test  section  at  G  =  1.19  x  10^  and  T^  ^  290°F, 
Accidental  burnout  and  time  limitations  precluded  going  further  with 
this  study. 
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Data  Reduction 

All  data  reduction  was  done  by  hand  calculation  with  the  aid  of  a 
desk  calculator.  This  included  computation  of  local  flow  and  heat  trans- 
fer parameters,  including  flow  rate,  heat  flux,  wall  superheat,  bulk  temp- 
erature or  enthalpy  as  well  as  the  dimensionless  groupings  Re,  Nu,  Pr, 
Nu/Pr^'^,  The  mass  velocity  (G)  was  determined  from  the  fluid  flow  rate 
and  tube  dimensions.  The  temperatures  were  converted  from  thermocouple 
emf s  with  NBS  standard  calibration.  The  heat  input,  was  obtained  from 
the  tube  voltage  and  current  (EI)  and  checked  by  the  resistance  and  cur- 
rent, a  thJ-rd  value  was  obtained  from  the  mass  flow,  specific  heat,  and 
temperature  rise.  Values  obtained  by  each  were  in  excellent  agreement 
and  differed  at  most  by  10^,  but  were  usually  within  $%   of  each  other. 
The  latter  method  was  of  no  value  when  the  exit  temperature  reached  sat- 
uration. The  average  heat  flux  was  calculated  using  the  several  values 
of  heat  input. 

Heat  transfer  results  were  based  on  measurements  of  the  outer  wall 
temperature,  A  correction  for  temperature  drop  through  the  tube  was  made 
to  first  obtain  the  inner  wall  temperature*  A  modified  form  of  the  Krieth 
and  Summerfield  solution  was  used  to  obtain  the  necessary  correction. 
This  solution  assumes  an  infinite  cylindrical  resistor,  with  heat  trans- 
fer at  the  inner  surface  and  an  adiabatic  outer  surface.  The  final  form 
of  the  equations  as  ivell  as  thermal- property  information  is  given  in 
Reference  (1)»  The  final  form  of  the  curve  (^^T^  versus  I  at  various  values 
of  T  J  is  shown  in  (Figure  3).  It  should  be  noted  that  at  high  values 
of  current,  the  wall  temperature  drop  is  extremely  large.  The  data  taken 
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in  the  last  run  (Table  IV  Series  N  Appendix  C )  points  out  the  effects 
of  this  on  the  results.  The  wall  temperature  drop  is  so  large  that  the 
inside  wall  temperature  is  actually  reduced  as  the  power  level  is  in- 
creased. Results  then  indicate  a  substantial  reduction  in  the  degree 
of  vxall  superheat.  This  leads  to  the  conclusion  that  the  equation  for 
wall  temperature  drop  may  not  produce  accurate  results  at  high  values 
of  current.  High  currents  were  reached  with  T.S.-2  as  well,  (Table  IV 
Series  F)  but  with  a  shorter  test  section  the  outside  wall  temperatures 
were  also  higher,  hence,  the  effect  was  not  as  pronounced.  At  each 
thermocouple  position,  calculations  were  made  to  obtain  the  local  heat- 
transfer  results.  The  fluid  bulk  temperatures  are  based  on  a  linear 
temperature  gradient  where  data  is  for  subcooled  bulk  fluid.  Where  net 
vapor  existed,  the  bulk  temperature  is  the  saturation  temperature  for  the 
particular  local  pressure.  The  local  pressure  was  obtained  graphically 
by  drawing  a  smooth  curve  through  the  observed  pressures  at  three  loca- 
tions. All  fluid  properties  are  based  on  these  local  boilk  temperatures. 
With  the  saturation  temperature  (Tg2.t)>  "t-^®  degree  of  subcooling  (A^g  v" 
■^  at  ~  '^b^  ^^^  ^^®  ""^■'"■'"  superheat  (ATg^^  »  T^  -  T   )  were  available, 
A  sample  set  of  calculations  is  shown  in  Appendix  B,  A  tabulated 
summary  of  the  data  obtained  is  given  in  Appendix  Co 
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RESULTS 
Critical  Heat  Flux 

Burnout  data  were  required  to  set  an  upper  limit  on  the  heat  flux 
that  could  safely  be  applied  to  the  heat  transfer  section.  This  data 
is  presented  graphically  in  Figure  U,  for  G  -  I.19  x  10^  Ib/hr  ft^ 
(V  c:^5ft/sec).  Also  plotted  are  data  taken  by  Wessel  (3)*  for  the  same 
tube  size  with  varying  fluid  velocity  and  a  slightly  larger  tube  with 

varying  L/D(V  ~  5ft/sec).  The  data  shows  an  increase  in  (q/A)   as 

cr 

exit  quality  increases. 
Heat  Transfer  Results 

The  non-boiling  results  of  this  study  are  presented  in  Figure  $, 
along  with  the  line  indicating  McAdaras  (U)  correlating  equation.  Data 
in  Figure  h   were  obtained  for  G  »  2.38,  1.19  and  0.59  x  10^  Ib/hr  ft^ 
(V:=:i(10,  $,  and  2,5  ft/sec).  The  span  of  Reynolds  number  also  results 
from  the  bulk  temperature  variation.  The  data  is  consistently  above 
that  predicted  by  McAdams  equation  and  can  be  represented  by: 

Nu  -  0.017  Re°*^^  Pr^'^  (2) 

Subcooled  nucleate  boiling  data  is  plotted  in  Figure  6.  Figure  6 
contains  data  at  the  same  flow  rates,  P  —  60  psia,  L/D  ■  37.i;»  and  a 
range  of  subcooling  from  about  0  to  175°F.  Three  other  curves  are 
plotted  to  aid  in  describing  the  results. 

First,  the  curve  of  (q/A)p-,  is  the  curve  of  heat  flux  versus  wall 
superheat  (T  -T   )  obtained  from  non- boiling  prediction.  This  curve 
is  based  on  G  «  1.19  x  10^  Ib/hr  ft^,  and  Tg^^  -  109-183°F,  where 


*  Numbers  denote  reference  listed  in  Bibliography 
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(q,/A)p^  -  ^(  V\^  "  ^^  V^sat  ^  "^sat  "  "^b^  ^^"^  ^  ^^  determined  from 
properties  at  bulk  temperature  corresponding  to  the  data  using  Eq  (2). 

Second,  the  curve  of  (q/A)^^  represents  the  locus  of  points  at  which 
surface  boiling  begins  and  where  the  heat  flux  departs  from  the  non- 
boiling  prediction.  Tnis  is  discussed  in  i^eference  (5)  and  is  given  by: 

/  0.023U 
(q/A),  -  15.60  p  ^-^^^  (T^  -  T^at)''  (3) 

The  curve  plotted  here  is  for  p  -  60  psia. 

The  third  curve  represents  the  mean  line  throuth  data  of  Wessel 
for  fully-developed  subcooled  nucleate  boiling  and  vias  determined  as 
indicated  in  Figure  (7). 

The  data  covers  a  range  of  flov7  conditions  in  the  transition  region 
from  the  inception  of  boiling  to  a  point  on  the  fully-developed  boil- 
ing curve  at  which  bulk  boiling  begins.  As  fully-developed  boiling  is 
reached,  the  data  approaches  an  asymptote  which  agrees  vri.th  livessel's 
data.  Figure  (6)  illustrates  that  in  the  transition  region  heat  trans- 
fer rates  are  dependent  on  ooth  flow  rate  and  subcooling,  and  as  the 
fully-developed  boiling  region  is  approached  the  dependence  on  flow 
rate  and  subcooling  disappears. 

Incipient  boiling  does  appear  to  begin  at  the  point  predicted  by 
the  intersection  of  the  incipient  boiling  curve  and  the  non-boiling 
prediction. 

Bulk  boiling  data  is  presented  in  Figures  8  and  9»  Figure  8  shows 
data  for  L/D  «  37. U,  at  the  two  lower  flow  rates  for  bulk  boiling 
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(i.e.,  at  points  where  the  bulk  fluid  was  at  the  saturation  tenpera- 
ture)o  The  data  represents  the  range  of  steam  quality  from  0  to  2k%m 
The  fully-developed  curve  is  again  plotted  for  comparison.   The  data 
of  Figure  8  to  appear  to  follow  the  fully-developed  boiling  curve  for 
the  range  of  quality  indicated. 

Figure  9  contains  data  obtained  with  T.S.  -2,  3,  and  3  modified. 
The  range  of  steam  quality  is  extended  to  almost  UO^,  and  data  is 
plotted  for  varying  L/D^s.   The  data  for  a  particular  L/D  can  still 
be  represented  by  a  form  of  Eq.   (1),  But,  as  quality  increases  (rep- 
resented by  increasing  L/D),  the  fully-developed  boiling  curves  are 
shifted  toward  the  lower  wall  superheat.  The  effect  of  quality  seems 
to  override  the  pressure  effect,  i.e.,  the  lower  pressure  results  in 
higher  wall  superheat. 

Similar  results  appear  in  Figure  10  for  data  with  subcooled 
boiling  conditions.  As  L/D  increases,  at  first,  wall  superheat  in- 
creases, and  then  decreases.   This  appears  to  be  independent  of  the 
degree  of  subcooling  within  the  range  covered  in  Figure  10  (  T    « 
30  -  150°F)„ 

In  an  attempt  to  determine  the  effect  of  steam  quality  on  heat 
transfer  coefficient.  Figures  11  through  Ik   show  the  experimentally  de- 
termined heat  transfer  coefficient  versus  steam  quality.  Figure  11  con- 
tains data  obtained  with  T.S.  -2  for  G  -  0.59  and  1,19  x  10  Ibj/hr  ft^. 
The  results  indicate  that  heat  transfer  coefficients  are  dependent  on 
heat  flux,  and  flow  rate  alone  in  the  range  of  steam  quality.  At 
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quantities  above  10-15^  the  data  appears  to  merge  for  each  flow  rate 
and  the  dependence  upon  heat  flux  disappears.  The  effect  of  flow  rate 
is  quite  pronounced  as  would  be  expected.  At  lower  flow  rates  the  effect 
of  quality  on  the  local  transfer  coefficient  becomes  apparent  at  lower 
values  of  quality. 

Figure  11  shows  the  results  of  T.S.  -2  and  3  at  G  »  I.I9  x  10^,  and 
Figure  12  shows  similar  results  for  G  -  0.^9  x  10^,  Data  obtained  from 
these  test  sections  appears  to  be  inconsistent.  Results  of  T.S.  -2  in- 
dicate an  increase  in  h  as  flow  proceeds  downstream.  If  the  data  ob- 
tained from  T.S.  -3  is  then  compared  with  this,  it  would  appear  that  h 
is  relatively  constant  between  37. U  L/D  and  57.9  L/D.  Bej^nd  that  point 
h  reaches  a  maximum  and  then  decreases.  Because  of  the  inconsistency  in 
the  data,  T.S,  -3  was  modified.  The  thermocouples  located  upstream  on 
T.S,  -3  (L/D  -  57.9  and  6U.1)  were  moved  to  positions  which  corresponded 
to  the  downstream  thermocouple  positions  of  T.S.  -2  (L/D  -  30.9  and  37, U), 

Data  taken  with  this  test  section  is  shown  in  Figure  13,  Also  plot- 
ted for  comparision  is  data  taken  with  T.S.  -2  and  3  at  G  -  1.19  x  10 
Ibs/hr  ft  ,  Data  for  each  test  section  was  taken  at  approximately  equal 
values  of  heat  flux  and  inlet  temperature,  in  so  far  as  data  was  avail- 
able from  the  other  test  sections.  Results  of  these  last  tests  compare 
favorably  with  that  obtained  from  test  sections  2  and  3»  The  absolute 
values  of  h  at  corresponding  L/D's  are  about  the  same.  The  only  differ- 
ence is  that  data  from  the  modified  test  section  shows  an  increase  in  h 
at  the  last  then:nocouple  position  instead  of  the  decrease  noted  with 
ToS.  -3 
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In  order  to  compare  the  bulk  boiling  data  with  previous  work  in 
this  area,  Figure  lU  shows  the  dimensionless  ratio  h/h^  versus  the 
Martinelli  parameter  (l/X^^),  where  the  Martinelli  parameter  is  defined 
by: 

•  1/X^^  .  (x/l-x)0-9  {p^  /pf'^     {M^/A4^f'^  (U) 

h  is  the  experimental  heat  transfer  coefficient  and  hj_is  the  liquid  heat 
transfer  coefficient  as  given  by  the  McAdams  correlation.  Also  plotted 
here  is  Chen's  F  function  and  his  prediction  of  h  at  parsuneters  approx- 
imating the  mean  values  for  the  data.  The  curve  lies  below  the  data  at 
the  lower  L/D's  of  T,S.  -2  ,  and  3  modified,  but  above  the  data  for  high 
L/D's  of  T.S.  -3  and  3  modified,  Chen's  correlation  is  presented  in  a 
later  section. 

Figure  15  shovjs  similar  data  at  G  «=  ,59  x  10  along  with  Chen's  F 
.function  and  his  predition  of  h  at  parameters  approximating  the  mean 
values  for  the  data.  Here  Chen's  F  factor  seems  to  divide  the  data  for 
low  L/D's  and  high  L/D's.  At  the  low  flow  rate  Chen's  prediction  of  h 
is  higher  than  the  experimental  values. 
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DISCUSSION  0?  RESULTS 


Critical  Heat  Flux 


In  the  course  of  obtaining  this  data  several  interesting  phenomena 
occurred.  Although  the  primary  intent  of  this  study  was  to  obtain  heat 
transfer  data,  it  is  considered  necessary  to  mention  these  items  briefly. 

First,  the  fact  that  (q/A)^^,  increased  as  much  as  ^0%   as  the  exit 

quality  increased  from  0  to  about  2^%.     Previous  studies  in  this  area 

id-th  subcooled  fluid  conditions  at  the  exit  have  indicated  that  (q/A) 

decreases  as  the  amount  of  exit  subcooling  decreases.  The  most  recent 

studies  show  that  as  zero  subcooling  is  approached  (q/A)   appears  to 

cr 

level  off  and  in  some  instances  increases  .  This  is  indicated  by  the 
data  of  Wessel  shown  in  Figure  U.  One  of  the  best  known  studies  of 
burnout  with  exit  quality  is  that  of  Lowdermilk  (6).  His  results  in- 
dicate that  (q/A)^,j.  is  at  a  maximum  when  exit  quality  is  approximately 
20^  depending  on  other  flow  parameters  as  well.   The  data  did  not  cover 
the  low  range  of  exit  quality.  Hence,  previous  work  indicated  that,  at 
a  constant  pressure  and  flow  rate, (q/A)   decreased  as  the  exit  subcool- 
ing  decreased.   Critical  heat  flux  reaches  a  relative  minimum,  increases 
in  the  low  range  of  exit  quality  to  a  relative  maximum,  and  then  falls 
off  sharply,   3ergles  (1)  also  shows  an  increase  in  critical  heat  flux 
with  an  increase  in  quality  up  to  12^,  Hence,  the  data  presented  here 
is  in  agreement  with  other  studies.  This  critical  heat  flux  is  not  as 
high  as  that  attained  with  subcooled  exit  conditions.  For  more  infor- 
mation in  this  area  the  reader  is  referred  to  recent  studies  by  Loosmore 
and  Skinner  (7)  on  critical  heat  flux  with  subcooled  boiling  and  by 
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Lopina  (8)  on  critical  heat  flux  :,dth  exit  quality  conditions. 

Second,  the  actual  location  of  burnout  of  the  test  section  was 
not  alvjajs  at  the  exit.  In  fact,  at  exit  quality  near  zero,  and  with 
exit  quality  about  2$%^   burnout  did  occur  at  the  exit.   Between  these 
limits,  the  location  of  burnout  was  considerably  upstream  of  the  exit. 
In  one  run  vrLth  exit  quality  approximately  13.6^  burnout  occurred  at 
18.7  L/D,  while  the  overall  tube  dimensions  in  all  runs  was  Itl.U  L/D, 
Calculations  indicate  that  the  bulk  fluid  condition  at  the  burnout 
location  was  slightly  subcooled.   This  indicates  that  in  the  low  range 
of  exit  qualities,  exit  conditions  may  not  be  the  correct  parameter  to 
consider  in  determinin?:  critical  heat  flux.  Lopina  (8)  and  VJaters,  etal, 
(9)  have  also  observed  upstream  burnout  and  discuss  this  in  more  detail. 
Waters'  study  shows  that  upstream  burnout  is  reproduceable  and  not  neces- 
sarily caused  by  instability. 

Third,  it  iira.s  observed  during  operation  that  as  the  exit  quality 
increased,  the  pressure  drop  from  the  test  section  exit  to  the  exit 
plenum  chamber  became  quite  large.  Hence,  1^Jhen  measuring  critical 
heat  flux  with  net  vapor  generation,  exit  pressure  must  be  determined 
at  the  test  section  exit  and  not  in  an  exit  plenum  chamber.  Results 
of  this  study  indicated  that  there  was  essentially  no  pressure  dif- 
ference between  these  two  locations  at  zero  quality,  but  at  17.8^ 
quality  the  pressure  difference  >ra.s  approximately  12.5  psia  at  G  ■ 
1.19  X  10°  Ibs/hr  ft^.  To  explain  this  it  is  necessary  to  give  a 
brief  discription  of  the  method  used  to  control  the  pressure. 

The  exit  pressure  level  was  controlled  by  means  of  a  plug  valve 
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dovmstream  of  the  exit.  At  low  values  of  heat  flux,  this  valve  was 
about  midway  between  the  fully  opened  and  closed  positions  in  order 
to  maintain  the  desired  pressure.  As  the  heat  flux  was  increased, 
and  boiling  conditions  were  reached,  the  exit  pressure  began  to  in- 
crease. This  valve  was  then  opened  to  reduce  the  pressure  to  the 
desired  level.  When  bulk  boiling  conditions  were  reached,  it  vras 
observed  that  pressure  control  was  no  longer  possible,  since  the 
valve  was  wide  open.  As  indicated  previously,  when  the  longer  test 
section  was  used,  pressure  control  was  lost  prior  to  obtaining  the 
maximum  heat  flux  for  that  series  of  runs.  This  then  illustrates  the 
choking  effect  that  causes  the  pressure  level  to  increase  in  the  test 
section* 

As  the  steam-water  mixture  leaves  the  test  section,  it  enters  a 
plenum  chamber  made  up  of  3/U  inch  brass  pipe.  Expansion  of  the  steam 
can  take  place  immediately.  The  pressure  drop  encountered  from  the 
test  section  exit  and  the  plenum  chamber  is  the  result  of  this  rapid 
expansion. 
Heat  Transfer  Results 

Non-boiling  data  Sxhown  in  Figure  5  plots  consistently  above  the 
correlation  of  McAdams,  The  divergence  appears  to  be  primarily  due  to 
the  radial  property  variation  which  is  not  properly  accounted  for  in 
the  correlation.  This  temperature  difference  effect  is  clearly  demon- 
strated in  Reference  (10).  In  any  case,  Eq.  (1)  is  not  intended  to  be 
a  generally  valid  correlation,  but  is  merely  indicated  to  represent  the 
data. 

Figure  6  shoxis  the  portion  of  the  boiling  curve  from  the  inception 
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of  boiling  to  the  fully-developed  boiling  curve  for  three  flow  rates 
and  two  ranges  of  subcoolings.  As  indicated  previously,  inception  of 
boiling  did  occur  at  the  predicted  value  of  wall  superheat  and  heat 
flux.   Beyond  this  ooint,  the  slope  of  the  curve  increases  until  fully- 
developed-boiling  is  reached.  Using  the  data  of  Vessel  to  establish 
the  fully-developed  boiling  curve  (Figure  7),  the  data  does  approach 
this  line  as  an  asymptote.  Kence,  the  subcooled  data  appears  to  be  con- 
sistent with  the  vrork  of  previous  studies. 

At  the  onset  of  this  study.  Professor  Bergles  pointed  out,  that  at 
high  heat  fluxes  with  low  subcooling,  the  data  appeared  to  be  shifted 
from  the  fully-developed  curve,  toward  lower  wall  superheats.  It  has 
already  been  mentioned, that  fully-developed  subcooled  boiling  data 
(Figure  10)  show  a  change  in  wall  superheat  at  different  positions  a- 
long  the  length  of  the  tube.  At  the  inlet,  subcooled  liquid  velocity 
can  be  specified  as  Vq_»  As  the  fluid  moves  down  the  heated  sections, 
surface  boiling  begins  and  non-equilibrium  voids  are  formed.  Therefore, 
the  volume  occupied  by  the  liquid  is  reduced.  In  order  that  continuity 
be  maintained,  the  liquid  velocity  must  increase.  Let  this  new  velocity 
be  V2»  As  boiling  becomes  more  vigorous  along  the  tube,  more  volume  is 
occupied  by  vapor  and  the  liquid  velocity  continues  to  increase.  This 
illustrates  how  the  liquid  velocity  may  increase  along  the  length  of  the 
tube.  ViJith  this  in  mind, and  noting  from  Figure  6  that  by  increasing 
fluid  velocity,  fully-developed  boiling  is  reached  at  higher  values  of 
heat  flux. 

If  we  consider  three  fluid  velocities  which  do  not  differ  by  much, 
and  determine  the  heat  flux  at  which  the  lowest  velocity  reaches  the 
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fully-developed  boiling;  curve,  it  is  clear  that  at  this  same  heat 
flux,  the  wall  superheat  will  be  reduced  at  the  two  higher  velocities. 
This  example  shows  how  the  wall  superheat  can  decrease  at  higher  L/D's 
\-ilth   subcooled  liquids.  This  also  indicates  that  the  local  heat  trans- 
fer coefficient  must  increase  in  this  range.  The  heat  transfer  co- 
efficient is  defined  as: 

h  -  (q/A)  /  (T^  -  T^)  (5) 

As  long  as  the  bulk  fluid  is  subcooled,  the  bulk  temperature  is  in- 
creasing  along  the  direction  of  flow.  If  (T^  -  Tg^^)  decreases  as 
indicated  previously  and  Tgg^^  decreases  with  pressure,  then  T^  must 
also  be  decreasing.  As  a  result,  (T^  -  T,  )  must  decrease  and  the  heat 
transfer  coefficient  increases  due  to  the  presence  of  void  fractions. 
This  explains  the  trend  tovxard  lower  wall  superheat  in  data  at  high 
L/D's  as  shown  in  Figure  10. 

Bulk  boiling  data  (Figure  9)  also  shox^s  a  reduction  in  wall  super- 
heat with  increasing  L/D.  Bennett,  etal,  (11)  presents  a  discussion  on 
the  effect  of  net  vapor  generation,  similar  to  the  discussion  above  on 
the  effects  of  void  fraction  on  fluid  velocity.  This  explains  the  lower 
wall  superheats  obtained  at  the  exit  of  the  test  section.  In  the  region 
of  bulk  boiling  where  the  bulk  temperature  is  the  saturation  temperature, 
a  decrease  in  wall  superheat  means  an  increase  in  h  for  a  specified  heat 
flux.  Therefore,  wherever  net  vapor  exists,  there  is  an  increase  in  the 
local  heat  transfer  coefficient. 

Now,  consider  a  section  along  the  tube  where  the  bulk  fluid  is  slight- 
ly subcooled.  As  flow  proceeds  down  the  tube  the  bulk  fluid  approaches 
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saturation  conditions.  The  rate  at  which  bubbles  nucleate  at  the 
surface  decreases  since  the  fluid  velocity  is  greater.  At  the  same 
time,  the  bubbles,  vjhich  have  Just  departed  the  surface,  do  not 
collapse  and  the  volume  fraction  of  vapor  remains  nearly  constant. 
Then,  the  fluid  velocity  will  also  remain  constant  and  the  degree 
of  wall  superheat  predicted  by  the  boiling  curve  will  not  change. 
Therefore,  in  this  region  near  zero  quality,  the  heat  transfer  co- 
efficient is  relatively  constant  as  well.  Farther  along  the  heated 
section,  net  vapor  generation  continues  by  means  of  surface  boiling 
and/or  evaporation  at  the  liquid  vapor  interface.  At  this  point,  the 
fluid  velocity  begins  to  increase  again  and  continues  to  increase 
until  the  test  section  exit  is  reached.  This  appears  to  explain  the 
data  of  Figure  9,  as  well  as  that  of  Figure  11, 

Figure  12  with  data  of  T„S.  -2  and  3  at  G  -  1,19  x  10^  Ibs/hr  ft^ 
shows  the  discrepancy  in  the  data  obtained  from  these  two  test  sections. 
Data  obtained  vjith  T.S,  -3  shows  that,  as  x  increases,  h  reaches  a  max- 
imum and  then  decreases.  This  could  not  be  explained  adequately.  Con- 
sidering data  for  an  approximate  equal  heat  flux  for  the  two  test 
sections,  there  also  appears  to  be  a  discontinuity  in  the  data  between 
L/B  -=  39.1;  and  57 ,9.  The  value  of  h  increases  as  L/D  increases  to  39. U, 
but  apparently  there  is  no  further  increase  in  h  to  L/D  »  57.9.  Beyond 
this  point,  h  reaches  its  maximum  value  and  decreases.  Results  for  G  - 
5.9  X  10^  Ibs/hr  ft^  shown  in  Figure  12  are  similar. 

It  was  expected  that  data  obtained  with  T.S.  -3  modified  would  re- 
solve this  discrepancy.  Figure  13  shows  data  for  one  flow  rate  and  one 
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value  of  heat  flux  for  T.S.  -2,  3,  and  3  nodified.   The  results  of  this 
show  that  the  absolute  magnitudes  of  h  remained  about  the  same.  The  only 
difference  which  did  appear  consistently  was  that  at  high  L/D»s,  h  in- 
creased rather  than  decreased.   The  heat  transfer  coefficient  then  show- 
ed and  increase  over  the  first  half  of  the  T.S. ,  then  remained  a  rela- 
tively constant  value,  and  finally  increased  again  near  the  exit.   One 
possible  explanation  is;  at  the  lower  L/D's  flow  conditions  had  juit 
reached  the  point  at  which  the  void  fraction  is  nearly  constant.  Over 
a  short  range  beyond  this  point,  h  is  constant  with  x.   If  this  were 
true,  the  results  as  indicated  in  Figure  11,  are  misleading  and  quality 
does  not  significantly  affect  h  until  it  reaches  a  value  of  about  20!^, 

Another  method  of  evaluating  the  data  would  be  to  draw  the  mean 
line  through  the  data  from  both  test  section.  It  is  uncertain  v^ich 
data,  if  any,  is  actually  incorrect.  If  this  were  done  in  Figure  12, 
it  would  indicate  the  effect  of  quality  becomes  apparent  at  about  1$%, 
but  its  overall  effect  would  not  be  as  great  as  indicated  in  Figure  11, 

Figures  1$   and  16  show  the  same  data  on  another  set  of  coordinates, 
h/hj-  versus  l/X  ,  Dengler  and  Addoms  (12)  were  the  first  to  use  these 
coordinates  in  evaluating  this  type  of  data.  Here  h-  is  defined  as  the 
conventional  liquid  heat  transfer  coefficient  determined  by  McAdams* 
correlation.  They  proposed  a  correlation  for  two  phase  heat  tri=Jisfcr 
coefficients  in  terms  of  l/X   (the  Martinelli  parameter)  and  an  F 
factor  to  account  for  the  degree  of  nucleate  boiling  present.  As  the 
fluid  progresses  into  the  annular  flow  regime,  with  continuiously  in- 
creasing steam  quality,  the  increased  velocity  of  the  two  phase  mixture 
induced  by  the  vaporization  process  suppresses  the  nucleate  boiling  proc- 
ess. Beyond  this  point  the  heat  transfer  coefficient  becomes  governed 
by  the  forced  convection  process. 
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V7ith  this  in  mind  Chen  (13)  proposed  the  following  correlation 
for  the  tvo  phase  heat  transfer  coefficient: 

^  "  hmic  -^  ^lac  (6) 

where  hj^^  is  the  irdcroconvective  portion  due  to  the  boiling  mech- 
anism and  hence  dependent  on  liquid  and  bubble  parameters: 
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W  '   •  00122    ^ 1 ! 
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and  h^g^^   is  the  macroconvective  given  by: 

Again  hj^  is  given  by  McAdams'  equation.  F  is  a  function  of  the 
Martinelli  parameter  and  is  defined  as  the  ratio  of  the  tv:o  phase 
Reynolds  number  to  the  liquid  Reynolds  number,  S  is  called  the 
suppression  factor  and  approaches  unity  at  zero  flow  rate,  and  zero 
at  infinite  flow  rate.  Gouse  (li;)  suggests  that  this  suppression 
factor  should  go  to  zero  much  sooner  and  thus  acknowledge  the  possi- 
bility of  total  suppression  of  nucleation  at  finite  Reynolds  numbero 
Both  S  and  F  are  presented  graphically  in  Reference  (13) o 

Chen  restricts  the  applicability  of  this  correlation  to  systems 
with  annular  or  annular  mist  flow  in  vertical  geometries.  This 
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correlation  successfully  predicts  the  data  of  numerous  independent 
studies  to  -within  plus  or  minus  15^,  and  seems  to  be  one  of  the  best 
available.  Gouse  (lii)  also  states  that  there  is  no  reason  why  it 
should  not  apply  to  annular  flow  in  horizontal  tubes. 

Data  obtained  with  thermocouples  at  the  top  and  bottom  of  the 
tube  wall,  at  the  same  location  along  the  length,  produced  no  indica- 
tion of  stratification.  Thermocouples  indicated  that  the  temperature 
difference  obtained  in  all  cases  was  less  than  the  accuracy  of  the 
measurement.  Also,  the  greater  temperature  fluctuated  between  the  top 
and  bottom.  VJith  this  in  mind  Chen's  correlation  is  considered  a 
suitable  means  of  checking  the  data. 

As  noted  earlier  Chen's  prediction  seems  to  plot  betx-reen  the  data 
for  high  and  low  L/D's  indicating  that  the  difference  can  be  considered 
as  scatter  in  the  data.  It  must  be  noted,  that  scatter  of  such  a  con- 
sistent nature  is  unusual* 
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CQNCLUSIONS 

Critical  Keat  Flux; 

The  critical  heat  flux  increased  as  much  as  S0%   as  exit  quality 
was  increased  from  0  to  2^^, 

Hoat  Transfer 

1.  The  local  heat  transfer  coefficient,  obtained  at  subcooled 
boiling  conditions,  does  appear  to  increase  due  to  the  pre- 
sence of  voids  formed  during  the  boiling  process, 

2.  As  bulk  fluid  approaches  saturation  conditions  the  local  heat 
transfer  coefficient  is  constant  over  a  low  range  of  steam 
quality.  The  exiient  of  this  range  increases  as  the  flow  rate 
increases. 

3.  For  the  flow  rates  considered  in  this  investigation,  the  local 
heat  transfer  coefficient  shows  a  significant  increase  as  qual- 
ity is  increased  beyond  1S%„ 

U.  Chen's  correlation  appears  to  be  £iirly  accurate  in  the  low 
range  of  steam  quality,  but  it's  accuracy  improves  as  quality 
increases  and  nucleate  boiling  is  suppressed. 
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RECOKMSNDATIONS 

Critical  Heat  Flux 

1.  tJhen  bulk  fluid  conditions  at  the  exit  contain  vapor,  the 
exit  conditions  must  be  determined  by  the  pressure  at  the 
test  section  exit  rather  than  the  pressure  in  an  exit  plenum 
chamber. 

2,  Since  upstream  burnout  is  possible  -with  exit  quality  conditions, 
and  stable  flow,  it  may  be  of  interest  to  determine  if  some  para- 
meter, other  than  exit  quality  might  give  a  more  meaningful  in- 
dication of  the  results. 

Heat  Transfer 

1.  Additional  experiments  should  be  conducted  with  test  sections 
similar  to  T<,S,  -3  and  instrumentation  covering  a  greater  por- 
tion of  the  test  section.  This  is  suggested  to  see  if  there 
is  actually  a  region  along  the  test  section. where  the  local 
heat  transfer  coefficient  is  constant. 

2.  Additional  study  of  heat  transfer  coefficient  in  the  low  qual- 
ity range  is  still  required,  in  order  to  develope  a  more  accu- 
rate prediction  in  the  range  where  heat  transfer  rates  are  de- 
pendent on  the  combined  effects  of  boiling  and  convcctive  heat 
transfer. 
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APPENDIX  A 

SXPEiS-ISNTAL  APPAPJiTUS 

Flow  Loop  and  Poi-rer  Supply 

The  flow  loop  used  in  the  experimental  program  was  one  already  in 
operation  at  the  M.I.T.  Heat  Transfer  Laboratory.  This  facility  is 
shown  schematically  in  Figure  1.  It  was  a  closed  loop  system  in  which 
all  components  are  of  corrosion  resistant  materials.  The  system  con- 
tained the  main  circulating  pump,  an  accumulator,  the  test  section  line 
\idth  parallel  flow  meters,  preheater,  test  section  and  its  instrumenta- 
tion, the  bypass  line,  and  a  heat  exchanger  utilizing  city  water.  Aux- 
iliary equipment  included  a  fill  pump,  supply  tank,  degassing  tank,  and 
a  continuous  demineralizer.  The  test  section  power  vras  supplied  by  an 
a.c.  motor  -  d.c.  generator  set.  A  discription  of  the  major  components 
follows.  More  detailed  information  is  given  in  Reference  (1)» 

The  main  circulating  pump  was  a  two-stage,  turbine-type  pump  driven 
by  a  3hp  induction  motor.  It  orovidad  a  head  of  2^0  psi  at  3.6  gpm,  A 
bladder  type  accumulator,  pressurized  in.th  nitrogen,  located  at  the  pump 
exit  served  to  darap  out  pressure  fluctuations.   The  ball  valve  in  the  by- 
pass line  was  used  to  control  the  by-pass  flow  rate  and  pressure,  which 
in  turn  determined  the  pressure  in  the  test-section  line.  The  test  sec- 
tion line  contained  a  Fischer-Porter  flow  meter  and  a  set  of  preheaters 
in  series.  There  are  four  5kvj  preheaters  plus  an  additional  5kw  emersion 
unit  fabricated  in  a  3  inch  copner  tube  and  provided  vrith  0  to  100^  con- 
trol by   a  Power stat  auto-transformer.  This  arrangement  allovred  varying 
preheat  from  0  to  25kw,  one  being  variable  and  the  other  allowing  incre- 
ments of  5kw  each. 
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The  test  section  po>'er  was  supplied  by  rriotor-£;enerators.  The 
generators  were  driven  by  iiUO  volt,  3  phase  synchronous  motors.   Two 
36kw  d.c.  generators,  each  nominally  rated  at  12  volts  and  3000  amps 
were  connected  in  series.   The  generators  .'ere  provided  vdth  water 
cooled  shunts  in  parallel  v.dth  the  test  section  which  allowed  them  to 
be  run  open  circuited  at  the  test  section  for  starting,  or  after  burn- 
out. 
Instrumentation 

Instrumentation  was  available  for  reading  pressure  levels,  differ- 
ential pressure,  temperature,  test  section  flow  rate  and  test  section 
voltage  and  current.  Pressure  levels  were  read  on  bourdon-tube  gages 
as  shoT'm  in  Figure  1,  At  the  test  section  inlet  and  exit  v;ere  200  and 
100  psi  test  gages  xirith  an  accuracy  of  plus  or  minus  l/Ii^  of  full  scale. 

Test  section  flow  rate  was  determined  from  the  larger  flowrater  meter 
which  had  interchangable  tubes  and  floats.  The  ranges  covered  are  20- 
160,  35-UOO,  and  lUO-1360  Ib/iir.  All  units  were  calibrated  as  installed 
in  the  system. 

The  power  supplied  to  the  test  section  was  obtained  from  the  test 
section  voltage  and  current.  The  voltage  was  read  directly  on  a  Vfeston 
multiple-range  d.c.  voltmeter  ;ri.th  a  specified  accuracy  of  plus  or  minus 
\/2%,     The  current  was  determined  from  a  standard  shunt  v/ith  a  calibra- 
tion of  60.11   amp/mv. 

All  temperatures  were  measured  with  copper-constantan  thermocouples 
made  from  30  gage  duplex  T-ri.re,  The  fluid  bulk  temperature  at  the  inlet 


i 
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•and  exit  of  the  test  section  was  measured  by  inserting  thermocouples 
directlj'  into  the  fluid  using  conax  fittings  with  lava  sealants. 
Methods  of  measuring  v/all  temperature  are  deecribed  later.  The  out- 
put voltages  of  the  thermocouples  and  the  shunt  v/ere  displayed  on  a 
continuous  recorder.  The  recorder  was  a  Brown,  single  channel  in- 
strument having  ranges  of  0-6,  5-11,  10-16,  15-21,  and  20-26  mv. 

The  manometer  system  used  to  measure  differential  pressure  con- 
sists of  two  Meriam  60  inch  U-tube  manometers,  manifolds  and  valves, 
and  connecting  lines  of  rubber  hose.  The  system  is  designed  to  read 
a  maximum  of  10  pressure  differences  with  either  one  or  tv;o  reference 
points  on  either  manometer.   One  manometer  was  filled  ivith  mercury 
which  gave  a  maximum  range  of "approximately  25  psi.  The  other  mano- 
meter contained  an  oil  with  a  specific  gravity  of  2.00  which  allowed 
a  maximum  pressure  difference  of  approximately  2  psi.  A  typical  test 
section  with  three  pressure  taps  is  sho-i'in  in  Figure  2.  The  pressure 
drop  from  the  reference  "A"  located  at  the  inlet,  to  points  2,  and  3 
could  be  read  on  either  the  mercury  or  oil  manometers  depending  on  the 
magnitude  of  the  pressure  drop.  The  two  pressure  gages  were  connected 
to  the  inlet  and  exit  to  give  the  pressure  level  and  a  check  on  the 
pressure  drop.  Vents  vrere  included  so  that  all  lines  could  be  purged 
of  air. 
Measurement  of  Tube-Wall  Temperature 

The  most  difficult  problem  in  any  heat  transfer  experiment  is  the 
measurement  of  the  temperature  of  the  heat  transfer  surface.  The  basic 
method  of  measuring  tube  wall  temperature  is  described  in  Reference  (l)o 
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The  method  used  here  v/as  as  xollo-.vs: 

Heated  shields  constructed  of  standard  aluminum  tubir^^  vrith 
threadlike  grooves  cut  into  the  outer  wall,  were  wrapped  i-dth  in- 
sulated Chromel  A-26  wire.  Spacers  fabricated  from  aluminum  served 
to  hold  the  shield  in  place. 

Power  was  suppled  in  the  Chromel  heating  coil  by  a  110  VAC  supply 
which  is  stepped  doxm  to  about  25  volts  by  a  Stancor  filament  trans- 
former. The  power  to  each  shield  is  controlled  by  a  Povierstat.   This 
arrangement  permits  line  control  of  shield  heating. 

A  layer  of  Scotch  Electrical  tape  Number  2?  was  wrapped  around 
the  tube  at  each  thermocouple  location.  The  measuring  junction  was 
then  bound  radially  around  the  tube  Td.th  another  layer  of  tape.  Thermo- 
couples were  also  taped  to  the  inside  of  the  shield.  Sufficient  length 
of  all  thermocouples  were  left  inside  the  shield  to  avoid  conduction 
errors.  The  space  between  the  shield  and  tube  was  filled  with  asbestos 
fiber.  Precautions  were  taken  to  insure  the  shield  was  electrically 
insulated  from  the  test  section, 

Povjer  to  the  shields  is  adjusted  imtil  the  shield  thermocouples 
and  the  tube  wall  thermocouples  within  the  shield  section  indicate 
approximately  equal  WiF» 
Test  Section 

The  test  section  was  constructed  from  30h   stainless  steel  tubing 
with  an  inside  diameter  of  0.2U25  inch.  Brass  bushings  were  cut  from 
3/U  inch  round  stock.  Holes  viere  drilled  for  the  test  section.  The 
bushings  were  then  sweated  to  the  tube  v;ith  silver  solder.  Holes  0.099 
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inch  were  then  drilled  in  the  bushings  for  pressure  taps.  These 
holes  did  not  penetrate  the  test  section.  A  number  75  drill  (0.021 
inch)  was  used  to  drill  the  pressure  tap  hole.   Beryllium  copper  pres- 
sure tap  tubes  (2irj7i  ID)  were  then  set  in  the  bushings  and  soft  soldered 
in  place.   The  inside  of  this  oube  was  filled  v:ith  milk  of  magnesia  to 
prevent  the  tap  being  plugged.  V/ith  this  completed,  the  milk  of  magnesia 
could  then  be  removed. 

An  additional  pressure  tap  was  placed  near  the  center  of  the  test 
section.  Here  a  short  stainless  steel  tube  was  first  filled  with  milk 
of  magnesia  and  then  silver  brazed  to  the  tt:st  section.  Excess  braze 
was  then  filed  away  to  give  a  very  small  fillet  at  the  test  section. 
The  milk  of  magnesia  was  removed  and  the  pressure  tap  hole  drilled  (0.021 
inch)  through  the  test  section.  Any  burrs  inside  the  tube  were  removed 
with  fine  emery  cloth  and  pieces  of  steel  wool  pushed  through  the  tube, 

Thermocouples  were  then  attached  to  the  test-section  wall  as  des- 
cribed previously.  Thermocouples  were  located  at  four  positions  along 
the  length.  At  two  locations  the  thermocouples  were  placed  perpendicu- 
lar to  the  tube  axis,  and  the  lead  wares  wrapped  around  the  tube  and 
taped  in  place.  At  the  other  two  locations,  thermocouples  were  again 
placed  perpendicular  to  the  tube  axis,  but  at  each  of  these  positions 
two  thermocouples  were  used,  one  on  the  top  of  the  tube  and  the  other 
at  the  bottom.  This  set  up  v/as  intended  to  detect  stratification. 

The  guard  heaters  were  then  assembled  onto  the  test  section.  After 
the  first  guard  heater  was  in  place  a  section  of  0.099  inch  OD  tubing 
was  placed  over  the  center  pressure  tap  and  soft  soldered  to  the  stain- 
less steel  tube  to  extend  this  pressure  tap,  A  diagram  of  this  test 
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section  is  shovm  in  Figure  2, 

ToS,  -1  and  2  Here  identical  in  construction.  -Accidental  burnout 
early  in  this  investigation  is  the  reason  for  this.  The  location  of 
pressure  taps  and  wall  therraocouples  is  indicated  belovo 

T.S,  1  and  2 

Gaining  length   6  in.  Distance  to  thermocouples 

Heated  length   10  in.  No,  1        lii.U  L/D 

Center  pressure  tap  at  27.3  L/D        2        23.7 

3       30.9 
h  39.1 

T.S.  3 

Calming  length   6  in.  1  57.9  L/D 

Heated  length   20  in.  2  6U.1 

Center  pressure  tap  at  U$,6   L/D        3  72. u 

U  78.5 


T.S.  3  modified 
Calming  length   6  in. 
Heated  length   20  in. 
Center  pressure  tap  at  hS.S   L/D        3        72. U 

i;        78.5 


1  30.9  L/D 

2  39.U 
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app:lKdix  b 

SAI'^LE  C.hLCULATION 

1.  Calculat-ion  of  heat  input  to  the  test  section; 

(a)  from  the  current  and  voltage  measurements 
q^  »  3.ul3  EI  -  3.al3  (60.17)  E  Eg 

«  205.36  E  Eg  3tu/hr 

S  and  E  were  measured  values, 
s 

(b)  from  current  and  resistance 


I. 


q,   »  3.U3  l"-     ?t  «  1.236  X  10^  £^  R 
D  S 

R   «  p   L/a 
•  e     c 

For  T.S.  -1  and  2   L  =  10/l2  ft, 

^c  "=  (Tl/ii)  (ID^  -  0D^)/lJ4a  -  2.115x10^  ft^ 

R  -  Oo9U  X  10  f^ 
^Q   was  determined  from  Reference  (1)  as  a  function  of  the 
tube  temperature.  The  tube  temperature  is  given  by: 

^  =  ^ow  -  (AT-,,/2) 
AT  1-ras  determined  from  Figure  2  and  is  a  function  of  T 
lo  T  ^  was  measured  by  means  of  thermocouples  converted  with 
K.B.S.  tables, 

(c)  from  mass  flow  rate  and  bulk  temperature  change 

%     =  '^  %     ^'^exit  -  ^i) 

"^exit  ^^^  '^'    ^^   ''ell  a-3  all  other  temperatures  was  obtained 

as  indicated  above  for  T^^ . 

ow 

A  sample  set  of  data  and  calculations  are  shown  on  the  following  pages. 
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2ata  for  Run  A- 8 


S               Sg  Pj,  APi-2       A?i~3         P3  T.                  T  y.. 

V               rw  psig  in-oil       in-oil       psig  mv  ^F        mv         ^F 

12.65  18.62  i;6.o  0,9            8.9          lili.6  0.52  56      3.68      I89 

inv             nv  nv            F           ny             Of  nv             °F 

11.36  11.30  11.37  161      11.29      ii55.7  ll.UO  I462 


•«->\- 


mv  rav  mv  F  nv  °F  r.iv  '^F 

ll.iiS       11„1;3      11.66    U70.7     11.51    i;56.7  ll.ii3        1;63 

q^     »     12.65   (ie.62)    (205.36)     »  U8122     3tu/hr 

q^     -     lo236  X  10^  (18.62)2     1.119  x  10"^  «  17750     3tu/hr 

I     =     60,17    (18.62)     =     1120  amps  Sc  Tq^.,,  -  hS^  °F 

then  Figure  2  gives   A'^^  »  115  °F 

T.  =  I165  -  57.5  -  li.07.5  °F  and   ja^  =  2.61^  (I'Q-^)  ohm  ft 
R  =  l.n9  X  10""'^  ohms 
q^  =        36U   (iSy  -  ^6)      =     >t8)4l2     Btu/hr 

q  "   ^^a  ■"  Qb  "^  ^^c^  /^  "  ^^""'^^^  Btu/hr 
2.  Heat  transfer  surface  area  A  «  1T(ID)  L  «  0o0529  ft^ 

3-     Average  heat  fLjx     q/A  =  4B093/  .0529  -  9.159  x  10^     Btu/hr  ft'^ 

U,  The  change  in  enthalpy  through  the  test  section  is 

Ah'   •=  o/v  »  135  3tu/lb^ 

5.  The  change  in  bullc  temperature  from  inlet  to  any  position  along  the 


T,;_,  Tt^,  TQ-and  Tvj  denote  gur.rd  heater  temperatures. 

The  a  and  b  indicate  temperatures  at  the  top  and  bottom  of  the  tube 
respectively. 


tube  is  given  by: 

%  '   (Texit  -  Ti)  x'/  L 

6,  The  pressure  at  any  point  along  the  tube  \<us  obtained  graphically 
as  indicated  earlier, 

7,  Entering  the  stean  tables  T..dth  the  local  pressure,  the  saturation 
temperature,  the  saturated  liquid  enthalpy  and  the  latent  heat  of 
vaporization  vere  obtained, 

8,  All  fluid  properties  at  each  thcnr.o couple  position  vere  obtained 
at  the  local  bulk  temperature.  These  include  Cp,  k, />,>t(,  and  P  , 

9,  At  each  position  the  follovri-ng  values  were  obtained  as  indicated 
above . 


L/D   x»/L     T„^,    AT^    T     P  ,   T 

sa1 

X 


^"^^  ^'    '-J      X'    sat 


a  li6l,0 
lU.ii      .35  115      3U5.0      60.2      292.9 

b  1;58.7 

23.7      .575       1462,0      115      3U7.0      59.9      292.6 

30.9      .75         ii70.7      115      355.7      59.6      292.3 


a  1:65.7 
37.  U   .906    "^''"    115   3a9o3  ^"^.h       292.1 
b  l!63.0 


■-i-     The  temperatures  at  a  and  b  were  averaged  to  give  the  outside  wall 
temperature  at  these  locations. 
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L/2  AT>.  T,  AT     .  AT     ^  fT     -  T  ^ 

lii.U        i;6.3  102.3  52.1  191.6  2U2.7 

23.7         76  132  S$.k  160.6  215.0 

30.9        99  155  63  oU  lli7.3  200.0 

37»U  119o9  165..9  57.2"  126.2  l63.ii 

10.   The  locil  value  of  h  Viis  not  calculated  for  subcoolcd  boiling  data, 
but  was  determined  for  non-boiling  and  bulk  boiling  data  as  showi: 
h  »  (c/a)  /  (T^  -  T,^) 
llo  Nu  and  Re  were  determined  for  non-boiling  data  as  well  as  the  ratio 
Nu/Pr^°^  . 

Nu  =  hD/k     =  (h/k)     (0.0202) 

Re  =  GD/^(  «  2.U0i;  X  1(^  /ji{  ©     w  =  361;  Ibj^  /  hr 

Pr  =  G, 


P 
Pr  is  tabulated  as  a  function  of  temperature  in  Reference  (15)  with 

the  other  properties  for  watero 

12,  Bulk  boiling  data  included  some  additional  calculations  as  indicated; 

(a)  local  enthalpy  h«:^,   «  h»^  +  A^'  ,i 

where  A^-\r   =  (^^'^^^t  "  ^'i^  ""'/" 

NOTE:-  v-Jhen  this  value  «xceeded  the  value  of  saturated  liquid 

enthalpy,  net  vapor  was  present.  The  local  quality  was 

calculated  using  the  local  enthalpy,  saturated  liquid 

enthalpy  and  latent  heat  of  vaporization. 

(b)  Evaluation  of  the  Martinelli  parameter  required  values  of  yi\ 
and  JO   for  steam.  yMy  '""^^  obtained  from  Reference  (l6)  and 
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f^^   was  obtained  from  the  steam  tables.  The  Martinelli 
parameter  was  then  calculated  for  each  data  point  at  which 
net  vapor  was  present.   The  Martinelli  parameter  is  defined 
as: 

(c)  To  obtain  the  ratio  h/hj^j  the  liquid  heat  transfer  coefficient 
(hi  )  was  calculated  at  each  of  these  data  points  with  McAdams 
equation; 

h^  -  0.023  (k/D)  ReO-^  Pr°*^ 
This  equation  vjas  used  because  it  provided  a  means  of  compar- 
ing the  data  vriith  Chen's  \cork  in  this  field.  His  value  of  h^^^ 
is  equal  to  hr ,  as  determined  above,  multiplied  by  his  F  factor. 


A  summary  of  the  data  is  contained  in  Appendix  C, 
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APPENDIX  C 
SUT-!!"1PJIY  OF  !jATA 


Table  I 


Burnout  Data 


Table  II     Non-boiling  Data 


Table  III     Subcooled  Nucleate  Boiling  Data 


Table  IV      Bulk  Boiling  Data 


Table  V      Range  of  Variables 


Page 
UO 

Ul 
69 


TABL5  I 
BURNOUT  DATA 

D  -  0.2u25  inch    L  -  10  inches    L/D  «  Ul.i;     V  -  5  ft/sec 


RUN 

^exit 

psia 

(o/A)^^^x  10-^ 
Btu/hr  ft^ 

(h'   .^-  h' 
^  exit 

3tu/lb 

sat^ 

(h'n  -  h'   ,  ) 
^  BO     sat' 

Btu/lb 

1 

30 

1.325 

20 

^ 

2 

^6 

l.i;2 

65 

U7 

3 

52 

2.135 

2hl 

•* 

U 

UO 

1.18 

52 

-30.1; 

5 

Ii2 

1.568 

109 

1.1 

6 

US 

1.638 

160 

20.1 

T.S.2 

72 

1.95 

296.It 

■» 

-X-  Burnout  occured  at  the  exit 
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TABLS  V 
RANGE  OF  V;j^IAHLES 


-1 

G  X  lO*^ 

q/A  X  10"^ 

°F 

iVhr  rt2 

Btu/hr  rt^ 

U3-57 

1.19 

.013-. 95U 

a,b 

1U9-153 

lol9 

,092-. 918 

a,b 

150-168 

1.19 

.125-o979 

a,b 

I8ii-187 

1.19 

.085-1.33 

a,b 

212-228 

1.19 

.306-1.68 

b,c,r 

2h9-289 

1.19 

.177-1.79 

C    t. 

5U-95 

2.38 

.319-1.35 

a,b 

152-156 

2.38 

.22  -1.007 

a,b 

55-63 

.59 

.236-.175 

a,b 

158-166 

.59 

.159-.5285 

b 

271-285 

.59 

.202-..625 

c 

276-293 

1.19 

.3185-. 878 

c 

285-291 

.59 

.2016-.  701; 

c 

281-290 

1.19 

.5165-1.169 

c 

■^   denotes  types  of  data  obtained 
a  -  non-boiling 

b  -  subcooled  nucleate  boiling 
•  c  -  bulk  boiling 
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